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I refer the reader to Ref. [Has a supplement to this brief account. 

1. FACTORIZATION OF THE INCLUSIVE QUARKONIUM 
PRODUCTION CROSS SECTION 

In heavy-quarkonium production in hard-scattering processes, two large momentum 
scales appear: the heavy-quark mass m and the typical momentum transfer in the 
hard scattering, which I will denote generically by px- One would like to separate the 
perturbative physics at these large momentum scales from the physics at smaller mo- 
mentum scales that is associated with nonperturbative heavy-quarkonium bound- 
state dynamics. It has been conjectured^ that, for the inclusive quarkonium pro- 
duction cross section at pT 3> m, one can achieve such a separation and that one 
can write the cross section in the following factorized form: 

a(if)=^F„(0|Of|0). (1) 

n 

The Fn are "short-distance coefficients." They are essentially the process-dependent 
partonic hard-scattering cross sections convolved with the parton distributions. The 
partonic hard-scattering cross sections depend only on the large scales m and pr, 
and they have an expansion in powers of ag. The quantities (0|O^(A)|0) are long- 
distance matrix elements (LDMEs) that are formulated in terms of the effective field 
theory nonrelativistic QCD (NRQCD). They give the probability for a heavy QQ 
pair with a certain set of quantum numbers to evolve into a heavy quarkonium H. 
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The LDMEs have a known scahng with v, the heavy-quark velocity in the quarko- 
nium rest frame. Hence, the factorization formula in Eq. ([1} is a double expansion 
in powers of as and v. In practice, the sum in Eq. ([T]) is truncated at a finite order in 
V. The LDMEs are believed to be universal (process independent) — an assumption 
that gives the factorization formula much of its predictive power. 

A key feature of NRQCD factorization is that quarkonium production can occur 
through color-octet, as well as color-singlet, QQ states. If one drops all of the color- 
octet contributions and retains only the leading color-singlet contribution, then one 
obtains the color-singlet model (CSM). The CSM is inconsistent for production 
processes beyond leading order in v, for example, production through P-wave QQ 
channels, because it leads to uncanceled infrared divergences. 



1.1. Status of a Proof of Factorization 

A proof of the NRQCD factorization formula in Eq. ([T]) is complicated because glu- 
ons can dress the basic production process in ways that apparently violate factoriza- 
tion, for example, by connecting one incoming hadron to another or by connecting 
an incoming hadron to the quarkonium. A proof of factorization would involve a 
demonstration that diagrams in each order in can be re-organized so that (1) 
all soft singularities cancel or can be absorbed into NRQCD matrix elements and 
(2) all coUinear singularities and spectator interactions can be absorbed into parton 
distributions. 

Nayak, Qiu, and Sterman pointed out that, in order to make the color-octet 
NRQCD LDMEs gauge invariant, one must be modify them from the original 
definitiorp] by including eikonal lines (path-ordered exponentials of path integrals 
of the gaug e field) that run from the points at which the QQ pair is created to 
infinity.'^^'^ These eikonal lines are essential at two-loop order in order to allow cer- 
tain soft contributions to be absorbed into the LDMEs. However, they do not affect 
existing phenomenology, which is at tree order or one-loop order in the case of the 
color-octet channels. A key difficulty in proving factorization to all orders in 
is the treatment of gluons that have momenta of order m in the quarkonium rest 
frame. Such gluons can communicate with the QQ pair through the exchange of 
soft gluons, which lead to infrared divergences. One could absorb these infrared 
divergences into the LDMEs by treating the gluon with momenta of order m as an 
eikonal line in each LDME. However, such a treatment would preserve the univer- 
sality of the LDMEs only if they are independent of the direction of the eikonal 
line (the direction of the gluon). In an explicit calculation at two- loop order, the 
dependence on the direction of the eikonal line cancels. However, it is not known 
if this cancellation generalizes to all orders. An all-orders proof of factorization is 
essential because non-factorizing soft-gluon contributions, for which is not small, 
could be important numerically at any order. 

Nayak, Qiu, Sterman have also pointed out that, if an additional heavy quark 
is approximately co-moving with the QQ pair that evolves into the quarkonium. 
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then there can be soft color exchanges between the heavy quark and the QQ pairP 
This process does not fit into the NRQCD factorization picture because it involves 
LDMEs that contain additional heavy quarks beyond the QQ pair that evolves 
into the quarkonium. The process is nonperturbative and, therefore, its rate cannot 
be calculated reliably. However, the process can be identified experimentally. Its 
signature is an excess of heavy-flavor mesons in a cone around the quarkonium of 
angular size of order mv divided by the quarkonium momentum. 



2. Quarkonium Production in Hadron-Hadron Collisions 

It has been known for many years that the cross sections differential in px for the 
production of the J/ip, '0(25'), Xcj, and ^ states at the Tevatron can be fit well by 
the predictions of NRQCD factorization at leading order (LO) in as. The existing 
fits have included the contributions from the LDMEs of leading order in v in the 
color-singlet and color-octet channels, and, in the S'-wave case, the two independent 
color-octet LDMEs of the first subleading order in v. The color-singlet LDMEs can 
be fixed by decay rates, but the color-octet LDMEs are treated as free parameters. 
While NRQCD factorization passes the minimal test of fitting the Tevatron cross 
sections, it is important to make more stringent tests by using the fitted values 
of the LDMEs to predict quarkonium production rates in other processes and to 
predict values for other observables, such as quarkonium polarizations. 

At LO in as, NRQCD factorization predicts that prompt J/tjj and tjj{2S) polar- 
izations at the Tevatron should be substantially transverse for values of the quarko- 
nium pt that are greater than rricP This prediction is compatible, within uncer- 
tainties, with the CDF Run I measurement, but is incompatible with the CDF 

o 

Run II measurement. The Run I and Run II results are also mutually inconsistent, 
but the CDF collaboration endorses the Run II results. The CDF Run II data for 
the polarizatiorP are also incompatible with the LO NRQCD prediction. In 

the case of T polarization, the and Dcff^ measurements are incompatible 

with each other and with the LO NRQCD prediction.^IH The discrepancies between 
theory and experiment with regard to polarization present a significant challenge 
to our understanding of the mechanisms of quarkonium production. Comparisons 
of predictions for the polarization of the J/'ijj with experimental measurements are 
greatly complicated by the presence of feeddown contributions from the ip(2S) and 
Xcj states. Measurements the polarization of the J/4> in direct production would 
be of considerable help understanding the production mechanisms. 

In the case of quarkonium production in the color-singlet channel, complete cor- 
rections at next-to-leading order (NLO) in as and real-gluon corrections at next-to- 
next-to- leading order in as (NNLO*) have been calculated. These calculations have 
revealed the surprising result that, at large pr, higher-order contributions can be en- 
hanced by more than an order of magnitude relative to lower-order cont ributio ns be- 
cause the higher-order contributions fall off less rapidly as px increases .^^^HSIj^ ^t^q 
case of prompt J/ ijj production at the NNLO* color-singlet contribution^ 
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lies well below the data, suggesting that most of the cross section comes from a 
color-octet contribution. In the case of prompt T production at CDF|IS]the NNLO* 
color-singlet contributiorPS! could explain the data by itself. However, the large theo- 
retical uncertainties still allow the pres enc e of a substantial color-octet contribution. 
In the case of J/^ production at STARTS and at PHENIX,tlSl the NLO and NNLO* 
color-singlet corrections are large,^!^ and the upper limit of the large uncertainty 
band for the NNLO* prediction is close to the STAR data at large pr- 

NLO calculations of the ^Sq and ^5*1 color-octet contributions to S-wave quarko- 
nium production at the Tevatron and the LHC show that the corrections to J/ip 
and ipi'^S) production are less than 25%, while the corrections to T production 
are less than 40%.!^ Recently, two groups have completed the first NLO calcula- 
tion s that include all of the color-octet channels through the first subleading order 
in ^ | 22 | 23 |rp]^g numerical results of the two groups for the short-distance coefficients 
are in agreement and show a very large negative K factor in the ^Pj channels. In 
Ref. [231 the LDMEs were extracted in a fit to both the CDF Run II datiSI and 
the HI photoproduction data,^2lli25i and the resulting LDMEs are not qualitatively 
different in size from those that were obtained in L O f its to the Tevatron data. In 
Ref. [22j fits were made only to the CDF Run II data,!^ and two linear combinations 
of LDMEs were determined. The fitting procedures that were used in Refs. [22land 
[23l differ in a number of respects, beyond the inclusion/exclusion of the HI data, and 
the LDMEs that were obtained differ substantially, in one case by an order of mag- 
nitude. Clearly, it will be necessary to examine the validity of the assumptions that 
are implicit in the fitting procedures before any conclusions can be drawn about 
the sizes of the LDMEs. In Refs. [22] and [23l the extracted values of the LDMEs 
were use d to make predictions for the prompt J/ip cross section differential in pT at 
CMS,^^ and good agreement was obtained. In Ref. [23l a prediction was also made 
for the J ftp cross section at PHENIX,!^ and good agreement was again obtained. 
The fact that good agreement with the CDF and CMS data was obtained using 
very different values of the LDMEs suggests that the hadroproduction cross section 
differential in pT is rather insensitive to the details of the production mechanism. 

NLO and NNLO* calculations of the color-singlet contribution to the J/ip and 
T polarizations reveal th at they change from transverse to longitudinal once one 
goes beyond Lo! ^'^ ! ^^ ^^ NLO corrections to the color-octet 5*- wave contributions 
to the J/ip polarization produce only small effects.^^ Given the large sizes of the 
NLO corrections to the rates from color-octet ^Pj channels, it is clearly necessary to 
compute the contributions to the polarization from these channels before a definitive 
comparison can be made between NRQCD factorization predictions and experiment. 

The large corrections that occur at NLO and NNLO* cast some doubt on the con- 
vergence of the perturbation series. As I have mentioned, the large corrections arise, 
at least in part, because the NLO and NNLO* contributions in certain QQ channels 
fall more slowly with pT than do the lower-order contributions. The px distribution 
can fall no more slowly than This dependence is achieved in the color-singlet 

channel at NNLO* and in the color-octet channels at NLO. Hence, no further kine- 
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matic enhancements are expected beyond these orders. However, the theoretical 
uncertainties are large in these orders because the kinematically enhanced contri- 
butions are, in effect, computed at the Born level. Kang, Qiu, and Sterman have 
suggested that this difficulty might be overcome by using fragmentation-function 
methods to reorganize the perturbative calculation according to the pT dependence 

on 

of the contributions. Such a reorganization might make it possible to compute 
more accurately the contributions that are most important numerically. 

3. Quarkonium production in ep collisions 

Until recently, it had been believed that NLO color-singlet contributions account 
quite well for the J/'0 photoproduction cross sections tha t have been measured 
at HERA, leaving little room for a color-octet contribution.l2IIS21 However, a new 
calculation of the NLO color-singlet contribution,"^"^ while confirming the analytic 
results of the previous calculations, reveals that a more reasonable choice of the 
renormalization/factorization scale yields a much smaller numerical value for the 
color-singlet contribution, thereby opening the possibility that there is an appre- 
ciable color-octet contribution. This possibility was substantiated by a complete 

no OA 

NLO analysis, including color-octet contributions, ^ in which it was found that 
colo r-octet contributions are necessary in order to obtain agreement with the HI 
data.l23E51 While the NLO cross section differential in px is in good agreement with 
the HI data, the NLO cross section differential in the J/tp energy fraction z devi- 
ates from the data at low and high z. However, these deviations are understood as 
arising from the presence of uncalculated resolved-photoproduction contributions at 
small z and the failures of convergence of the NRQCD velocity expansion and the 
perturbation expansion near z = 1, both of which require a resummation in order 
to obtain reliable predictions. 

The color-singlet contribution to the J /iJ; polarization in photoproduction is 
strongly affected by QCD corrections ,1^2^^ changing from largely transverse at LO 
to largely longitudinal at NL O. Com parisons of the NLO color-singlet contribu- 
tion with the HI and Zeus datePEll show that the color-singlet contribution alone 
cannot explain the observed polarization. A complete NLO calculation, including 
color-octet contributions, would be necessary in order to make a meaningful com- 
parison between the predictions of NRQCD factorization and the data. 

An LO calculation of J/ip production in deep-inelastic scattering at HERA is 
generally in good agreement with the data.l^ However, given the importance of 
NLO corrections in other quarkonium production processes, it would be prudent to 
wait for the calculation of the NLO corrections to deep-inelastic scattering before 
drawing any conclusions. 
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